Natural Resolution of the Proton Size Puzzle 
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We show that ofT-mass-shell effects arising from the internal structure of the proton provide a new 
proton polarization mechanism in the Lamb shift, proportional to the lepton mass to the fourth 
power. This effect is capable of resolving the current puzzle regarding the difference in the proton 
radius extracted from muonic compared with electronic hydrogen experiments. These off-mass-shell 
effects could be probed in several other experiments. 
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The recent, extremely precise extraction of the proton 
radius [l[ from the measured energy difference between 
the ^Py^^ and 25*^^2^ states of muonic hydrogen (H) 
has created considerable interest. Their analysis yields a 
proton radius that is smaller than the CODATA value 
(extracted mainly from electronic H) by about 4% or 5.0 
standard deviations. This implies [l[ that either the Ryd- 
berg constant has to be shifted by 4.9 standard deviations 
or that the QED calculations for hydrogen are insuffi- 
cient. Since the Rydberg constant is extremely well mea- 
sured, and the QED calculations seem to be very exten- 
sive and highly accurate, the muonic H finding presents 
a significant puzzle to the entire physics community. 

Our analysis is motivated by the fact that muonic hy- 
drogen is far smaller than electronic hydrogen and there- 
fore more sensitive to corrections arising from hadron 
structure. In particular, we consider the lowest order cor- 
rection associated with off-shell behaviour at the photon- 
nucleon vertex, showing that it can very naturally ac- 
count for the difference reported by Pohl et al. Since 
at the present state of development of hadronic physics 
it is not possible to provide a precise value for this cor- 
rection, our result may be viewed as a phenomenological 
study of the sensitivity of muonic hydrogen to important 
aspects of proton structure. It should spur further study 
of processes which could be sensitive to off-shell changes 
in proton structure. In alternate language, the explana- 
tion which we present may be viewed as a new contribu- 
tion from proton polarization that is not constrained by 
dispersion relations but which can be studied in systems 
other than the hydrogen atom. 

We begin with a brief discussion of the relevant phe- 
nomenology. Pohl et al. show that the energy difference 
between the 2P^^^ and 25^^^^^ states, AE is given by 

AE = 209.9779(49) - 5.2262r^ + 0.0347r3 meV, (1) 

where rp is given in units of fm. Each of the three coef- 
ficients is obtained from extensive theoretical work 0- 
Q, typically confirmed by several groups. Studies of the 
relevant atomic structure calculations and corresponding 
efforts to improve those have revealed no variations large 
enough to significantly affect the above equation Q . 



Using this equation, we see that the difference between 
the Pohl and CODATA values of the proton radius would 
be entirely removed by an increase of the first term on 
the rhs of Eq. (1) by just 0.31 meV=3.1 x 10'^° MeV, but 
an effect of even half that much would be large enough to 
dissipate the puzzle. Finding a new effect of about that 
value resolves the puzzle provided that the correspond- 
ing effect in electronic H is no more than a few parts in 
a million (the current difference between theory and ex- 
periment 3). An effect that gives a contribution to AE 

of the form ct^-^ (with m the lepton mass and M the 
proton mass) could therefore resolve the proton radius 
puzzle and cause no disagreement in electronic H. 

The search to find such an effect has attracted consid- 
erable interest. New physics beyond the Standard Model 
must satisfy a variety of low-energy constraints and so 
far no explanation of the proton radius puzzle has been 
found that satisfies these constraints pi]|- [l4{ . Attention 
has been paid to the third term of Eq. ([T|) jl5l| , with the 
result that its current uncertainties are far too small to 
resolve the proton radius puzzle [H, [T^l • 

We therefore seek an explanation based on the fact 
that the proton is not an elementary Dirac particle, and 
that many features of its interactions are still unknown. 
In particular, consider the electromagnetic vertex func- 
tion which must depend on all of the relevant invariants. 
For a proton of initial four-momentum p, the most gen- 
eral expression must include a term, dependent on the 
proton virtuality, that is proportional to — M"^ and/or 
p ■ 7Ar — M, where the subscript N denotes acting on a 
nucleon, and M is the nucleon mass. Such terms have 
been discussed for a very long time in atomic [1, and 
nuclear physics (T^-jll]- They have been of special con- 
cern in relation to the difference between free and bound 
deep inelastic structure functions measured in the EMC 
effect [l^- 22 1, nucleon- nucleon scattering [2p a nd elec- 
tromagnetic interactions involvin g nu cleons |24l. [25j , no- 
tably quasi-elastic scattering [26j j32{. 

Many possible forms [l^, include the effects of pro- 
ton virtuality; we consider three that could be significant 
for the Lamb shift. We write the Dirac part of the vertex 
function for a proton of momentum p to absorb a photon 
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of momentum q = p' — p. as: 
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where three possible forms are displayed. Other terms of 
the vertex function needed to satisfy the WT identity do 
not contribute significantly to the Lamb shift and are not 
shown explicitly. The proton Dirac form factor, Fi{—q^) 
is empirically well represented as a dipole Fi{—q^) = (1 — 
g2/A2)-2, (A = 840 MeV) for the values of -q^ = Q'^ > 
of up to about 1 GeV^ needed here. F{—q^) is an off- 
shell form factor, and A+(p) = {p ■ 7jv + M)/(2M) is an 
operator that projects on the on-mass-shell proton state. 
We use Oa unless otherwise stated. 

We take the off-shell form factor F{—q^) to vanish at 
= 0. This means that the charge of the off-shell proton 
will be the same as the charge of a free proton, and is 
demanded by current conservation as expressed through 
the Ward-Takahashi identity [H, . We assume 
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This purely phenomenological form is simple and clearly 
not unique. The parameter b is expected to be of the 
order of the pion mass, because these longest range com- 
ponents of the nucleon are least bound and more suscep- 
tible to the external perturbations putting the nucleon 
off its mass shell. At large values of \q'^\, F has the same 
fall-off as Fi, if 5 = 0. We take A = A here. 

We briefly discuss the expected influence of using 
Eq. The ratio, R, of off-shell effects to on-shell ef- 
fects, R ~ (p-7iv-A/) ^2^^ < A2) is constrained by 
a variety of nuclear phenomena such as the EMC effect 
(10-15%), uncertainties in quasi-elastic electron-nuclear 
scattering [26j . and deviations from the Coulomb sum 
rule [13]. For a nucleon experiencing a 50 MeV central 
potential, {p ■ -fN - M)/M ~ 0.05, so \q^ /b"^ is of or- 
der 2. The nucleon wave functions of light- front quark- 
models [111 contain a propagator depending on M^. 
Thus the effect of nucleon virtuality is proportional to 
the derivative of the propagator with respect to Af , or of 
the order of the wave function divided by difference be- 
tween quark kinetic energy and M . This is about three 
times the average momentum of a quark 200 MeV/c) 
divided by the nucleon radius or roughly M/2. Thus 
R ^ {p ■ IN - M)2/M, and the natural value of \q^ jb^ 
is of order 2. 

The lowest order term in which the nucleon is suffi- 
ciently off-shell in a muonic atom for this correction to 
produce a significant effect is the two-photon exchange 
diagram of Fig. [T] and its crossed partner, including an 




FIG. 1: Direct two-photon exchange graph corresponding to 
the liitherto neglected term. The dashed line denotes the 
lepton; the sohd hne, the nucleon; the wavy lines photons; 
and the ellipse the off-shell nucleon. 



interference between one on-shell and one off-shell part 
of the vertex function. The change in the invariant am- 
plitude, M.OS-, due to using Eq. ^ along with O^, to be 
evaluated between fermion spinors, is given in the rest 
frame by 
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where the lepton momentum is I ~ (m, 0,0, 0), the vir- 
tual photon momentum is k and the nucleon momentum 
p = (Af, 0,0,0). The intermediate proton propagator 
is cancelled by the off-mass-shell terms of Eq. This 
graph can be thought of as involving a contact interaction 
and the amplitude in Eq. (j4]) as a new proton polariza- 
tion correction corresponding to a subtraction term in the 
dispersion relation for the two-photon exchange diagram 
that is not constrained by the cross section data [s^ ]. 
The resulting virtual-photon-proton Compton scattering 
amplitude, containing the operator 7^7 v corresponds to 
the T2 term of conventional notation js^], [l^. Eq. 
is gauge-invariant; not changed by adding a term of the 
form A:^ k'^ /k"^ to the photon propagator. 

Evaluation proceeds in a standard way by taking the 
sum over Dirac indices, performing the integral over k^ 
by contour rotation, k^ —ik^, and integrating over the 
angular variables. The matrix element M. is well approx- 
imated by a constant in momentum space, for momenta 
typical of a muonic atom, and the corresponding poten- 
tial V = iM has the form V{t) = V[)5{r) in coordinate 
space. This is the "scattering approximation" Then 
the relevant matrix elements have the form Vb 1^*25 (0)|^, 
where '^23 is the muonic hydrogen wave function of the 
state relevant to the experiment of Pohl et al. We use 
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(am^) /(Stt), with the lepton-proton re- 
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duccd mass, nir. The result 
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(where f3 = Am? / h?) shows a new contribution to the 
Lamb shift, proportional to ra^ and therefore negligible 
for electronic hydrogen. Using leads to a vanishing 
hyperfine HFS splitting because the operator 7^ is odd 
unless ^ = 0. 

We next seek values of the model parameters A,6, ^ 
of Eq. ([3]), chosen to reproduce the value of the needed 
energy shift of 0.31 meV with a value of A of order unity. 
Numerical evaluation, using ^ = 0, A = A, shows that 
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52 (79MeV)2 
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leads to 0.31 mcV. If ^ is changed substantially from 
to 1 the required value of Awould be increased by about 
10%. If our mechanism increases the muonic Lamb shift 
by 0.31 meV, the change in the electronic H Lamb shift 
for the 2S-state is about 9 Hz, significantly below the 
current uncertainty in both theory and experiment 0. 
Should some other effect account for part of the proton 
radius puzzle, the value of X/b^ would be decreased. We 
also caution that other systems in which one might aim 
to test this effect could show sensitivity to the value of ^ 
or A in Eq. 

The other operators appearing in Eq. ^ yield similar 
results when used to evaluate MoS- Using Ob, gives a 
term of the T2 form with a Lamb shift twice that of Oa, 
and also a HFS term that is about -1/6 of its Lamb shift, 
so the value of X/b"^ would be decreased by 3/5. The use 
of Oc, gives a term of the Ti form and the same Lamb 
shift as Oq, as well as a HFS term that is -1.7 times its 
Lamb shift. In this case, the value oiX/lP' would be about 
-3/2 times that of Eq. ©. The HFSs arc small enough 
to be well within current experimental and theoretical 
limits for electronic hydrogen. Thus each operator leads 
to a reasonable explanation of the proton radius puzzle. 

It is necessary to comment on the difference between 
our approach, which yields a relevant proton polarization 
effect, and others [36[, [13] which do not. The latter use 
a current-conserving representation of the virtual-photon 
proton scattering amplitude in terms of two unmeasur- 
able scalar functions, Ti_2- Dispersion relations are used 
to relate Ti^2 to their measured imaginary parts. How- 
ever, terms with intermediate nucleon states are treated 



by evaluating Fcynman diagrams. This allows the re- 
moval of an infrared divergence by subtracting the first 
iteration of the effective potential that appears in the 
wave function. But the Feynman diagrams involve off- 
shell nucleons, so that their evaluation for composite par- 
ticles must be ambiguous. For example, using two differ- 
ent forms of the on-shell electromagnetic vertex function, 
related by using the Gordon identity, leads to results that 
differ. This ambiguity in obtaining Ti^2 is removed in our 
approach by postulating Eq. ^ and evaluating its con- 
sequences. Note also that in order to evaluate the term 
involving Ti using a dispersion relation one must intro- 
duce a subtraction function, Ti{Q,q^). This is uncon- 
strained by prior data [3^ because the value of ctl/ot at 
J/ = cxD is not determined Pachucki [s^, in Eq.(31), 
assumes a form proportional to (see our Eq. ([3])) times 
the very small proton magnetic polarizability. However 
we are aware of no published derivation of this result. 

In conclusion, we have shown that a simple off-shell 
correction to the photon-proton vertex, which arises nat- 
urally in quantum field theory and is of natural size and 
consistent with gauge invariance, is capable of resolv- 
ing the discrepancy between the extraction of the proton 
charge radius from Lamb shift measurements in muonic 
and electronic hydrogen. Off-shell effects of the proton 
form factor were an explicit concern of both Zemach Q 
and Grotch & Yennie However, it is only with the re- 
markable improvement in experimental precision recently 
achieved [l[ that it has become of practical importance. 
Within the field of nuclear physics there is great inter- 
est in the role that the modification of nucleon structure 
in- medium may play in nuclear structure [sol . |40| . We 
stress that the effect postulated here can be investigated 
in Icpton-nucleus scattering via the binding effects of the 
nucleon, as well as by lepton-proton scattering in arenas 
where two photon (or 7, Z) effects are relevant. 
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